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Motivated by recently reported experimental phase diagrams, we study the effects of CoOe 
distortion on the electronic structure in Na2,Co02 • 1/H2O. We construct the multiband tight- 
binding model by employing the LDA result. Analyzing this model, we show the deformation 
of band dispersions and Fermi-surface topology as functions of Co02-layer thickness. Consid- 
ering these results together with previous theoretical ones, we propose a possible schematic 
phase diagram with three successive phases: the extended s-wave superconductivity (SC), the 
magnetic order, and the spin-triplet SC phases when the Co valence number s is -1-3.4. A phase 
diagram with only one phase of spin-triplet SC is also proposed for the s=-|-3.5 case. 
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Since the discovery of superconductivity (SC) in 
Naa;Co02 • ?/H20,^"^ several experiments have suggested 
a spin-fluctuation-mediated unconventional pairing. 
However, its SC properties have been controversial, prob- 
ably because the sample preparation needs delicate con- 
trol. Even its pairing symmetry has not been clarified yet. 
Several experiments have provided results that are incon- 
sistent with each other, such as those for upper critical 
field i?c2,^ heat capacity^"' and phase diagrams. ^'^ 
This makes the situation more complicated. In partic- 
ular, phase diagrams proposed in the early stage of re- 
search are contradictory. One group proposed a constant 
Tc as a function of Na content, -'^^ while another group 
proposed a dome-shaped behavior of Tc^^ 

It has recently been pointed out that the valence of 
Co ions in this material is not -1-3.65 as naively expected 
from the Na content x of ~0.35.^'*'^^ Instead, the ac- 
tual valence is -1-3.4-3.5 irrespective of x because of the 
presence of oxonium ions (HsO"*"). Here, the valence s is 
directly related to the number of t2g electrons per Co ion 
(ritjg) as nt2g = 9 — s. Thus, the tuning of x does not nec- 
essarily correspond to the carrier doping. This fact has 
not been recognized so far, and the contradiction of ex- 
perimentally obtained phase diagrams mentioned above 
is, in fact, attributed to this point. 

Elaborate phase diagrams have been experimentally 
obtained by precisely determining the Co valence s by 
Sakurai et al.^'^^'^"^ They proposed a phase diagram in 
the plane of temperature and x. They found that when 
one prepares samples whose s value is constant at -1-3.4, 
a SC phase (SCI), a magnetically ordered phase (MO) 
and another SC phase (SC2) successively appear as x in- 
creases. On the other hand, for samples with s =-f3.5, 
only one SC phase appears. Here, we can consider that 
X controls the trigonal CoOg distortion along the c-axis 
(in other words, the thickness of the C0O2 layers, which 
is constructed from edge-shared CoOg octahedra): A 
larger- a; sample has thicker C0O2 layers. Several groups 
have shown the relationship between Tc and the crys- 



tallographic c-axis parameter c/a and that between Tc 
and the thickness of C0O2 layers.^®' In addition, simi- 
lar phase diagrams have also been proposed recently by 
other groups.^' 

These phase diagrams indicate that the electronic 
structures in this material are quite sensitive to the thick- 
ness of C0O2 layers. Therefore, in order to clarify the 
properties and mechanism of the SC as well as the mag- 
netic properties, it is important to consider the effects 
of CoOg distortion more seriously. In fact, the present 
authors previously examined these effects by tuning the 
trigonal crystal field (CF) parameter in the multiorbital 
Hubbard model. ^'^ It was shown that the increasing trig- 
onal CoOg distortion pushes the e'g bands up, resulting in 
the formation of hole pockets near the K-points, and thus 
the ferromagnetic and triplet-pairing instabilities are en- 
hanced. However, in this study, the changes in hopping 
parameters due to the distortion were not taken into ac- 
count, and Sakurai's phase diagram with three successive 
phases (SCI, MO and SC2 phases) was not reproduced. 

In this letter, we study the effects of CoOg distortion 
on the electronic structure in Naj;Co02 • yH2 by con- 
structing a multiband tight-binding model including not 
only the fivefold Co 3d orbitals but also the threefold 
O 2p orbitals. (Note that our previous model contains 
only the threefold Co t2g orbitals. ^■^) Using this model, 
we first study microscopically the changes in both band 
dispersions and Fermi surface (FS) topology. We show 
that when s is fixed at -1-3.4, three types of FSs succes- 
sively appear as the CoOe distortion increases. We dis- 
cuss that this FS deformation can explain well Sakurai's 
phase diagrams for s =-f3.4. In particular, we predict 
that two SC phases in this phase diagram have differ- 
ent pairing states, one is a spin-singlet and the other is 
a spin-triplet. We also study the band-structure defor- 
mation and the phase diagram for the s =+3.5 case. In 
addition, the contradictory results of Hc2^ and heat ca- 
pacity measurements as well as those for the nature 
of magnetic correlations are also discussed. 
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Fig. 1. (a) Co t2g band dispersion calculated from multiband 
tight-binding model for <f>oCoO=^7-5° (dots), and that from LDA 
calculation for bilayer- hydrate system (solid lines). ■^^ Three hor- 
izontal dashed lines denote the Fermi levels for s =+3.4, +3.5 
and +3.67 cases, respectively, from upper to lower, (b) Trigonal 
CoOe distortion along the c-axis. Arrows indicate the directions 
of O-ion shifts, (c) Original Brillouin zone and folded Brillouin 
zone for expanded unit cell (see text). 



The multiband tight-binding Hamiltonian is given by 
H'^P = H„y. + y,s,dl^di, 
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where d^ and are the creation operators of the Co 3d 
electron and the O 2p electron, respectively. The indices 
i and j run over the Co and O sites, respectively. The 
indices 7 and I run over the fivefold 3d and threefold 2p 
orbitals, respectively. The first term represents CFs. The 
second and third terms express the level energies of the 
Co 3d and O 2p orbitals, respectively. The other three 
terms express the d-d, d-p and p-p electron transfers, re- 
spectively. 

Here, the term i?cry. consists of two components: 

Hcry. = ^egdl^^dima+ Dm'm"dl^i^dim"aj 

i,m,(T i,m\m" ,a 

(2) 

where m in the first term runs over the twofold Co eg 
orbitals, and m' and m" in the second term are indices 
of Co t2g orbitals. The first term represents the cubic CF 
generated by the octahedrally coordinated O ions, which 
contributes to the t2g-eg sphtting. The second term rep- 
resents the; trigonal CF due to the trigonal CoOe dis- 
tortion where Dm'm" = ^^D3d{Sm'm" — 1); which con- 
tributes to the aig-e'g splitting of threefold t2g degener- 
acy. 

The transfer integrals tdd, tdp and tpp are deter- 
mined by the Slater-Koster (SK) parameters:^'' {dda) = 
-0.395, {ddn) = 0.103 and {dd6) = 0.026 for the 
nearest-neighbor (NN) Co-Co bonds, {pda) = —2.060 
and (pdir) = 0.989 for the NN 0-Co bonds, ippa)i = 
0.623 and (pp7r)i = -0.164 for the NN 0-0 bonds, 
{ppcr)2 = 0.467 and {pptt)2 = —0.123 for the second- 



NN 0-0 bonds, ippa)3 = 0.327 and (^^77)3 = -0.078 
for the third-NN 0-0 bonds, and {ppa)4, = 0.101 and 
ippTr)i = -0.024 for the fourth- and fifth-NN 0-0 bonds. 
In addition, the charge-transfer energy — Sp is 2.20, 
and the CF parameters A^^ and Ansd are 0.55 and 0.0, 
respectively. These parameter values are deduced by fit- 
ting the LDA result for bilayer-hydrate material. ^'"^ The 
energy unit is eV. 

The C0O2 layer becomes thinner when the octahcdra 
are more significantly contracted along the c-axis. We 
simulate this trigonal distortion by displacing the O ions 
along the trigonal directions (see Fig. 1(b)). With this 
distortion, the OCoO bond angle ((^oCoo) is no longer 
90°. The angle is increased further from 90° when the dis- 
tortion is more significant. When we carry out the above 
fitting, 0oCoO is fixed at 97.5° along with the condition 
for which the LDA calculation was performed. 

In Fig. 1, we display the Co t2g band dispersions cal- 
ctilated using these parameters together with those ob- 
tained in the LDA calculation, which shows good coin- 
cidence. We take the expanded hexagonal unit cell with 
lattice vectors of length \/3a, which results in a folded 
Brillouin zone (BZ) with one-third of the volume of the 
original one. The F', K' and M' points are symmetry 
points defined on the folded BZ, and the relationship be- 
tween the original and folded BZs is shown in Fig. 1(c). 

From the relationship between {ppa)i and {ppa)2, we 
see that (ppa) and (ppn) for adjacent O ions are pro- 
portional to ^oo^^' ■*^here loo is the distance between 
two adjacent O ions. In addition, from the expression 
of the interatomic matrix element, we expect that {pda) 
and (pdn) are proportional to Iq^Q' where IcoO is the 
neighboring Co-0 distance. Using these relationships 
as well as the obtained SK parameters, we investigate 
the changes in band structure due to CoOe distortion by 
varying ^oCoO- 

Before discussing the results, we note that the trigo- 
nal CF parameter AdscI necessarily changes with vary- 
ing <^oCoO- Although the value monotonically increases 
with increasing 0ocoOj its estimation is quite difficult 
and contains ambiguity because we need to estimate the 
value of the local dielectric constant for the Coulomb in- 
teraction between electrons on the O and Co ions. Thus, 
we first discuss the effects of transfer-integral variation 
due to the distortion by fixing ADsd to be constant, and 
then we discuss the eff'ects of Adsci variation by assuming 
several ways of monotonic change. 

Figure 2(a) shows the Co t2g band dispersions for sev- 
eral (/>ocoO values. Here, Ansd is fixed to be constant at 
zero. The horizontal dashed line denotes the Fermi level 
for s =+3.4. As shown in this figure, both the aig band 
at the r point and the top of the e'g bands near the K 
points are located below the Fermi level with thick C0O2 
layers of (/)ocoO=96.5° and 97.0°. As the C0O2 becomes 
thinner, both are pushed up. Thus, both appear above 
the Fermi level when the C0O2 layers are thin (see the 
cases of (/)ocoO=98.0° and 98.5°). 

We show the deformation of FS topology for both 
s =+3.4 and s =+3.5 cases in Figs. 2(b) and 2(c), re- 
spectively. We see that in the case of s =+3.4, FS with 
two aig-band cylinders (FSl) appears for thick C0O2 lay- 
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Fig. 2. (a) Band dispersions calculated from multiband tight- 
binding model for several </)oCoO values. Deformation of Fermi 
surface for s=+3.4 (b) and that for s=+3.5 (c) with varying 
</>OCoO- 
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Fig. 4. Density of states (DOS) for case C. The left panel shows 
the total DOS, the partial DOS for the aig orbital and that 
for the e'g orbital from upper to lower. The vertical dashed line 
denotes the Fermi level for s=-|-3.4. The right panel magnifies 
the shaded area in the left panel. 
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Fig. 3. Phase diagram of Fermi-surface topology for s=+3.4 in 
plane of 0oCoO and Ansd- 

ers, while FS with one oig-band cyhnder and six e^-band 
pockets (FS2) appears for thin C0O2 layers. On the other 
hand, in the case of s =-f 3.5, FSl does not appear even 
for a thick C0O2 layer. Instead, for a wide range of 0oCoO 
values, FS2 appears, while another type of FS with two 
a\g cylinders and six pockets (FS3) is realized when 
the layer is very thick with 0ocoO=96.5°. 

Next, we discuss the effects of trigonal CF, i.e., Ansd- 
In Fig. 3, we show a phase diagram of FS topology in the 
plane of 0ocoO and Aosd for s=-|-3.4. We find four types 
of FSs: FSl, FS2, FS3 and FS4. Here, FS4 consists of only 
one a\g cylinder. In the actual material, we expect that 
Ansd will monotonically increase with 0ocoO- Here, we 
consider several ways of increasing Aosd, for example, 
as represented by lines for cases A, B, C and D in this 
figure. In the prior calculation, Ansd is constantly fixed 
at 0.0, which corresponds to the line denoted as case A. 
If Aosd increases along the line of case B, two types of 
FSs (FSl and FS2) successively appear with increasing 
(/)ocoO- On the other hand, we see that three types of 
FSs, FSl, FS3 and FS2, successively appear for case C. 
For case D, in particular, the successive deformation of 



these three FSs occurs even within a tiny </)oCoO change 
of 1°. 

This phase diagram suggests that FSl tends to be re- 
alized in a system with thick C0O2 layers, while a sys- 
tem with thin C0O2 layers tends to have FS2. In fact, 
the former is the FS which Kuroki et al. employed to 
propose a possible extended s-wave pairing,^'' while the 
latter is the one on which several theoretical proposals 
of spin-triplet p- and /-wave pairings are based. As 
for the spin fiuctuation, finite-momentum (q =/=0) fluctu- 
ation is predicted for FSl,'^'^ while nearly ferromagnetic 
fluctuation at zero momentum (q —0) is predicted for 
pg2. 28-30 rpj^g ^ fluctuation for FSl is enhanced by 
the electron-electron scattering between inner and outer 
aig FSs owing to the intraorbital Coulomb interaction,^'' 
while the q —0 fluctuation for FS2 is caused by the inter- 
band scattering between aig and e'g bands owing to the 
Hund's-rule couphng.^^ 

For a moderate thickness of the C0O2 layers such 
as (j)ocoO—Q7 -5° , the FS topology strongly depends on 
Ansd- If Aosd is small, as in case A, FS4 is realized. 
With such a FS, the system has no speciflc enhanced 
spin fluctuation, as discussed in previous work.^'^ On the 
other hand, when Ansd is moderate or large, FS2 or FS3 
appears as in cases B, C and D. With these FSs, we ex- 
pect enhanced spin fluctuations.^'^' In particular, we 
expect that the spin fluctuation is markedly enhanced 
for FS3 because both electron scattering between the in- 
ner and outer aig cylinders and scattering between aig 
cylinders and e'g pockets cooperatively contribute to the 
enhancement of spin fluctuation. 

In Fig. 4, we display the total density of states (DOS), 
the partial DOS for the Co aig orbital (aig DOS) and 
that for the Co e'g orbital {e'g DOS) for case C. The blue 
dashed line denotes the Fermi level for s=-|-3.4. This fig- 
ure shows the following three aspects: [1] For thick C0O2 
layers (such as 0ocoO=96.5°), the aig DOS is rather 
large while the e'g DOS is small. [2] For thin C0O2 lay- 
ers (such as 4>ocoO—Q8.5°), the e'g DOS is rather large 
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Fig. 5. Schematic pliasc diagram for tlic s =+3.4 case in plane 
of ijioCoO ^'iid temperature. The inset shows a schematic phase 
diagram for the s =+3.5 case. 



while the aig DOS is small. [3] Both aig and e'g DOSs are 

large for moderate thicknesses such as (/)ocoO=97.5° and 
98.0°. The situation in the moderate case is favorable for 
the enhancement of magnetic fluctuation due both to the 
Hund's-rule coupling between the aig and e'g bands and 
to the intraorbital Coulomb interaction in the aig band. 

Now, we discuss the experimentally obtained phase di- 
agrams for s=+3.4 and s=+3.5 (see Fig. 5) on the ba- 
sis of the above results as well as previously obtained 
ones.^^"^'^ For s=+3.4, the observed two SC phases are 
likely to be the singlet extended ,s-wave pairing with 
FSl^"^ and the triplet /- or p-wave pairing with FS2,^^"^° 
respectively. The inbetween MO phase is caused by the 
cooperative contributions from enhanced q ^ and q = 
spin fluctuations. This is owing both to the FS topology 
of FS3 (or FS2) and to the DOS structure in which both 
aig and e'g DOSs are rather large. On the other hand, 
the experimental phase diagram for s=-|-3.5 consists of 
only one SC phase. The pairing state in this SC phase 
is expected to be the triplet p- or /-wave one with FS2, 
from Fig. 2(c). Note that the c-axis parameter c/a, NQR 
frequency Vq and Na content x appearing in Fig. 5 are 
scaled with the thickness of C0O2 layers.®' ^^'^^ 

We speculate that NMR experiments, which recently 
observed the decreasing Knight shift below Tc, were 
performed for samples in the singlet-pairing phase. 
We expect that a constant behavior may be observed 
in a sample with triplet pairing having thinner C0O2 
layers. It is also worth pointing out that the confus- 
ing results of Hc2^ and heat capacity^"' measurements 
may be attributed to the difference in pairing state de- 
pending on the sample. In addition, the contradictions 
of NMR/NQR and neutron-scattering results for the 
magnetic-correlation type may also be due to this sample 
dependence.^"®' ^^'^^'^^ Thus, it is highly desired to per- 
form these experiments by determining where the sam- 
ple is located in the phase diagram. Angle-resolved pho- 
toemission spectroscopy for the bilayer-hydrate system 
is also desired to detect the predicted FS deformation, 
since few measurements have been reported to date.'^^ 

To summarize, we have constructed the multiband 
tight-binding model including the full degeneracies of Co 
3d and O 2p orbitals for NaxCo02 • yH20. Using this 



model, we have studied the efl'ects of CoOe distortion 
on band dispersion, FS topology and DOS. On the basis 
of this analysis, we have proposed possible interpreta- 
tions of experimental phase diagrams for both s =-|-3.4 
and s =-|-3.5 cases. We have predicted that the following 
three phases successively appear in the s =+3.4 phase di- 
agram with decreasing Co02-layer thickness: a SC phase 
with singlet extended s-wave pairing, a MO phase and 
another SC phase with triplet p- or /-wave pairing. On 
the other hand, only one SC phase with triplet pairing 
has been predicted for the s =+3.5 phase diagram. 
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